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Abstract 

In  this  paper,  we  used  correlative  matrix  as  a  feature  criterion  in  K-L 
transformation  for  feature  compression,  and  then  the  synthetic  matched  filter 
(SMF)  was  prepared  by  means  of  the  ^nthetio  disorimlnant  function  (SDP).  This 
method  can  effectively  redut .  the  size  of  feature  image.  By  this  filter,  the  shift,  scale 
and  rotation  invariant  optical  pattern  recognition  was  realized  with  a  high  signal 
noise  ratio  of  correlation  output. 

Key  words:  information  processing,  optical  pattern  recognition,  matched  niter, 
computer-generated  hologram,  K-L  transformation, 

I.  Introduction 

Since  Vender  Lugt  [1]  proposed  using  the  holographic 
recording  method  in  fabricating  the  well-known  complex-number 
m-atched-^tatial-f ilter  great  progress  has  been  obtained  in 
coherent  optical  pattern  recognition.  However,  the  com^plex- 
number  matched-spatial-filter  is  incapable  of  solving  the  problem 
of  multi-aspect  invariant  pattern  recognition.  Later  on, 

Casasent  et  al  proposed  using  the  synthetic  discriminant  function 
(SDF)  to  fabricate  a  matched  filter  [2];  thus,  the  invariant 
deformation  problem  is  solved  more  satisfactorily.  However, 
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there  are  still  some  imperfections  to  be  solved  [3] .  One  of  them 
is  the  reduction  of  the  signal-to-noise  ratio  with  increase  in 
the  num.ber  of  training  tests  [4].  In  this  paper,  the  correlative 
matrix  is  used  as  the  feature  criterion  and  the  number  of 
characteristic  patterns  for  compression  in  the  K-L  transform  [5] . 
Then  the  SDF  method  is  employed  to  synthesize  the  pulse  response 
function  of  the  matched  filter.  By  using  the  method  of  computer¬ 
generated  holograms  [5] ,  a  synthetic  matched  filter  is 
fabricated.  In  the  coherent  optical  pattern  processing  system, 
three-aspect  (translation,  rotation,  and  scaling)  invariant 
optical  pattern  recognition  is  accomplished. 

II.  Filter  Algorithms 


1.  The  SDF  method 

Various  deformations  of  the  target  pattern  are  composed  into 
a  function  sequence  {/•}(»»— 1,  .  Let  us  assume  that  the 

pulse  response  h(x,y)  of  the  matched  filter  is  the  linear  set 


ff. 


h(jD,  y)^ 


<1) 


if  the  inputted  pattern  is  a  certain  deformation  f^^  of  the  target 
pattern,  it  is  required  that  the  output  of  the  correlation  of 
h(x,y)  at  the  origin  of  coordinates  is  a  constant,  for  example, 

jr(0,  y)h{w, 

By  substituting  Eq.  (1)  into  Eq.  (2),  we  obtain 

1,  ♦‘■••“JJ/.C®.  y'/mC®;  dy, 


(2) 


(3) 


Define  {f„}  as  a  correlative  matrix 


iJ. 


r  fix 

fu 

...  fxB  ■ 

^  ... 

r** 

» 

• 

...  fis 

• 

• 
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(4) 

r,. 

...  fsjl  - 

for 

the 

matrix 

Eq.  (3) 

can  be  written 

as  follows: 
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U,  ,  (5) 

In  the  equation  U  is  a  row  vector  for  the  matrix  element  at  (1); 
row  vector  a  is  determined  by  the  following  equation 

a=R-^U^  (6) 


From  the  equation,  the  function  h{x,V)  can  be  determined  as 
satisfying  Eq.  (2) ,  thus  obtaining  a  matched  filter  with  the 
deformation-invariant  feature. 


2.  K-L  Transform 

Briefly  speaking,  by  using  a  characteristic  vector 
corresponding  to  the  characteristic  value  of  the  matrix,  thus 
constituting  an  orthogonal  transform  (of  the  transform  matrix) , 
which  is  called  the  K-L  transform  [5] .  If  for  a  certain 
M-dimensional  vector  all  characteristic  values  of  the  covariant 
matrix  are  arranged  according  to  a  decreasing  sequence 
(Xi>X3>—  >X,),  the  transform  matrix  constituted  by  the 

characteristic  vector  corresponding  to  the  first  K 
characteristic  values,  to  this  M-dimensional  vector  the  K-L 
transform  is  connected.  Then  the  M-dimensional  vector  is 
compressed  into  a  K-dimensional  vector.  From  the  geometric 
significance,  it  can  be  considered  as  a  new  M-dimensional  space 
transformed  from  a  vector  of  M-dimensional  space.  In  the  new 
space,  the  projections  onto  M—K(1<K<M)  axes  is  zero  or 
very  small,  therefore  it  can  be  neglected.  Thus,  a  M-dimensional 
problem  is  transformed  into  a  K-dimensional  problem. 


3.  Algorithm  of  synthetic  matched  filter 

To  apply  the  K-L  transform  to  a  pattern  processing  problem, 
the  first  problem  to  be  solved  is  how  to  select  the  initial 
vector  before  conducting  the  K-L  transform,  thus  compressing  the 
number  of  dimensions.  The  method  used  is  as  follows: 


assume  that  there  exist  M  training  sample  patterns;  in  each 
pattern  there  are  N=LxL  sampling  points,  then  we  can  have  a 
vector  sequence  O'-l,  .  There  are  n  dimensions  in 

each  vector.  From  this  M  N-dimensional  vector  we  can  constitute 
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a  matrix  A 


...  ai*>- 

ai*> 

• 

• 

• 

• 

• 

• 

• 

o'** 

-  aT- 

(7) 


There  are  nxm  orders  for  this  matrix.  Let  B-A^  and  we  obtain 

1—  [Bi,  Bi,  Bjf],  (8) 


ai*’ 

...  a?’  ' 

ai*> 

—  <4** 

• 

• 

• 

• 

• 

• 

... 

Here,  (*-1,2, n)  is  of  the  M-diraensional  vector. 

This  N  vector  is  selected  as  the  initial  vector. 

From  the  exterior  product  of  (5i-J5,)  /  we  obtain  a  matrix  C 

«  -  (&) 

Here  is  the  probability  mean  of  vector  B^.  Assume  in  M 
training  samples  that  the  emerging  probability  is  the  same;  by 
selecting  The  matrix  C  is  of  MxM  order.  By  solving  the 

intrinsic  value  X  and  the  intrinsic  vector  W  of  this  matrix,  and 
arranging  the  intrinsic  values  in  a  decreasing  sequence,  we 
obtain  and  (i— li  2,  .Sfj  **“li  2, 

Here,  with  the  exception  of  are  all  zero. 

Thus,  it  makes  the  method  ineffective  in  compressing  the  number 
of  characteristic  patterns  according  to  the  magnitude  of 
characteristic  values.  As  the  characteristic  vector 
corresponding  to  zero  characteristic  values  is  not  unique,  a  new 

criterion  should  be  sought  after. 

Let  us  use  the  correlative  matrix  as  the  new  criterion. 
Thus,  first  we  derive  from  Eq.  (9)  the  intrinsic  value  W  for  the 

transform.  v  • 

^CS)  I 

L  (10) 
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thus  we  obtain 


In  the  equation  (j-l,  2, —,if;  2,  •••,!»), 

a  matrix  B' 

rof’  *•* 

a;«)  •••  a2*' 


(II) 


-1'- 


■  o?*’ 

With  transposition  of  B'  and  letting  obtain 

r<i<«  -  4*’- 

flia)  _ji)  ...  MM) 

**  (12) 

•  •  • 

•  •  •  ' 

.aji)  a«*)  ...  af»i. 

Each  row  vector  of  A'  is  considered  as  a  new  characteristic 

pattern;  there  is  a  total  of  M  patterns.  Next,  solve  for  the 

correlative  matrix 


-rn 

ru 

rn 

• 

rai 

« 

• 

-fjfi 

• 

rjfj 

Ajj  Is 

-1.  2, 

(13) 


The  self -correlated  A-.  is  arranged  in  a  decreasing  sequence; 


K  values)  A!^': 

this  new  K  characteristic  pattern,  a  synthesis  is  conducted  by 
using  the  SDF  mentioned  previously;  we  can  obtain  a  pulse 
response  function  of  the  synthetic  matched  filter  with  high 
signal-to-noise  ratio  with  the  deformation-invariant  feature. 


To 


Since  the  characteristic  vector  corresponding  to  the  zero 
value  is  not  unique,  possibly  all  values  on  the  main  diagonal  of 
the  correlative  matrix  are  greater  than  the  threshold  value,  thus 
incapable’  of  effectively  accomplishing  the  characteristic 
compression.  At  this  stage,  a  special  solution  should  be  re¬ 
established,  then  solving  for  the  homogeneous  equation  set  in 
order  to  conduct  characteristic  compression  by  again  using  the 
characteristic  criterion  method  of  the  correlative  vector  until 
we  obtain  a  satisfactory  characteristic  compression. 
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III.  Experimental  Results  and  Discussion 

To  accomplish  the  three-aspect  (translation,  rotation, 
scaling)  invariant  optical  pattern  i:*ecognition  by  using  the 
synthetic  matched  filter,  first  calculate  the  filter  by  using  the 
algorithm  described  in  this  paper.  Thus,  by  selecting  letter  "E" 
as  the  target  pattern,  with  three  different  scale  ratios 
{0.5:1:1.5)  and  five  orientations  (within  the  range  between  zero 
degrees  and  60  degrees,  with  a  15-degree  spacing),  totalling  15 
patterns  as  the  initial  sample  for  the  K-L  transform, 
characteristic  compression  is  conducted  by  using  the 
characteristic  criterion  method  of  the  correlative  matrix.  As 
revealed  in  the  calculation  results,  in  this  problem  the  15 
characteristic  patterns  can  be  compressed  into  three 
characteristic  patterns.  Then  there  are  three  characteristic 
patterns  as  used  a  new  function  sequence;  by  using  the  SDF 
method,  synthesize  the  matched  filter. 

Before  conducting  the  optical  correlative  pattern 
recognition  experiment,  first  a  computer  was  used  in  simulation. 
The  letter  "E"  was  used  as  the  input  pattern  (refer  to  Fig.  1 
(a))  with  scale  ratio  of  1,  orientation  angle  of  0°,  as  well  as 
1.4  as  the  scale  ratio  and  30°  as  the  orientation  angle,  the 
synthetic  matched  filter  from  the  computation  was  used  as  the 
spatial  filter;  the  corresponding  three-dimensional  diagram  of 
input  power  spectrum  is  shown  in  Fig.  1(b).  For  convenient 
comparison,  the  SDF  method  is  also  used  to  synthesize  the  15 
initial  samples,  thus  obtaining  the  SDF  of  another  spatial 
filter;  the  three-dimensional  diagram  of  the  correlative  output 
power  spectrum  is  shown  in  Fig.  1(c).  From  Fig.  1,  we  can  see 
for  the  SDF  of  the  15  initial  samples,  the  top  terminals  of  the 
correlative  output  are  comparatively  flat;  thus,  it  is  difficult 
to  determine  the  position  of  maximum  value  as  the  background 
noise  is  very  high.  By  using  only  the  synthetic  matched  filter 
of  three  characteristic  patterns,  the  correlative  output  has  a 
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clear  sharp  peak;  the  background  is  relatively  low.  According  to 
the  following  equation,  calculate  the  signal-to-noise  ratio: 


8NB 


(®.  y)^dy 


$ 


a 


(14) 


In  the  equation,  g(x,y)  is  the  correlative  value  at  point  (x,y) 
of  the  output  plane;  S  is  the  area  of  the  entire  output  plane. 
Refer  to  Table  1  for  the  computational  results.  In  the  table. 

No.  1,  No.  2,  and  No.  3  indicate  the  input  diagram  with  scale 
ratios  of  0.5,  1,  and  1.5,  respectively.  We  can  see  that  within 
the  synthetic  range  of  the  synthetic  filter,  the  signal-to-noise 
ratio  of  the  synthetic  matched  filter  (SMF)  is  higher  by  about  20 
to  23  percent  of  the  SDF. 


IHg.  1  Computer  simulation  of  the  optical 
oorrelatiTe  calculation,  (a)  Input  samples, 
(b)  Output  results  with  SMF,  (o)  Output 
results  with  SDF. 


Fig.  2  -.Eiperiment  of  optical 


correlation:  (a)  input,  (b) 


output. 
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Finally,  by  using  the  method  of  computer-generated 
holograms,  the  calculated  SMF  is  made  into  a  Fourier-transf orm 
hologram  [6]  of  two  elements  Lohmann  III  type,  to  be  used  as  the 
spatial  filter  in  the  coherent  optical  pattern  processing  system. 
The  input  diagram  shown  in  Fig.  2(a)  is  placed  in  the  input 
plane,  then  the  output  results  obtained  on  the  output  plane  are 
shown  as  in  Fig.  2(b) .  From  Fig.  2  we  can  see  that  in  the 
optical  correlative  pattern  recognition  experiment,  if  target  "E" 
appeared  in  the  input  plane,  whatever  the  occurrence  of 
translation,  rotation,  or  scaling  variation,  only  the  scale  ratio 
and  orientation  angle  fall  within  the  synthetic  range  of  the  SMF, 
then  a  clear  correlative  bright  spot  appears  at  the  corresponding 
position  of  the  output  plane.  However,  for  other  letters  in  the 
input  plane,  there  are  no  clear  bright  spots  in  the  output  plane. 
This  indicates  that  the  three-aspect  (translation,  rotation,  and 
scaling)  invariant  optical  pattern  recognition  can  be  realized  by 
using  an  SMF. 

(1)  In  the  K-L  transform,  the  correlative  matrix  is  used  as 
the  characteristic  criterion,  thus  enabling  effective  compression 
to  be  carried  out  on  the  number  of  characteristic  patterns. 

(2)  The  characteristic  patterns  obtained  after 
characteristic  compression  are  used  as  a  new  function  sequence; 
the  fabricated  SMF  can  have  a  higher  signal-to-noise  ratio  than 
such  filters  fabricated  by  using  the  SDF  method. 

(3)  The  computer  simulation  and  optical  experimental  results 
prove  that  the  above-mentioned  SMF  has  a  three-aspect 
(translation,  rotation,  scaling)  invariance  in  the  correlative 
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output  within  the  range  of  the  synthetic  scale  and  orientation 
angle . 
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